Exact solution for the dynamical decoupling of a qubit with telegraph noise 
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We study the dissipative dynamics of a qubit that is afflicted by classical random telegraph noise 
and it is subject to dynamical decoupling. We derive exact formulas for the qubit dynamics at 
arbitrary working points in the limit of infinitely strong control pulses (bang-bang control) and we 
investigate in great detail the efficiency of the dynamical decoupling techniques both for Gaussian 
and non-Gaussian (slow) noise at qubit pure dephasing and at optimal point. We demonstrate that 
control sequences can be successfully implemented as diagnostic tools to infer spectral proprieties 
of a few fluctuators interacting with the qubit. The analysis is extended in order to include the 
effect of noise in the pulses and we give upper bounds on the noise levels that can be tolerated in 
the pulses while still achieving efficient dynamical decoupling performance. 

PACS numbers: 03.65.Yz, 03.67.Pp, 05.40.-a 



I. INTRODUCTION 



Small superconducting circuits containing Josephson 
junctions have received a great deal of attention re- 
cently as promising candidates for scalable quantum bits 
(qubits). Unfortunately, despite the last few years of 
remarkable experimental and theoretical breakthroughs 
00 000 neither of these circuits satisfies the re- 
quirement that the coherence times must be considerably 
longer than the gate operation times Thus a major 
challenge for the next years is to achieve control of the 
decoherence afflicting these devices. 

It is generally recognized that in superconducting 
qubits the main sources of noise are: the external cir- 
cuitry, the motion of trapped vortices in the supercon- 
ductors and the motion of background charges (BCs) in 
associated dielectrics and oxides. While the noise orig- 
inating from the external circuitry can be reduced and 
practically eliminated by improving the circuit design 
and the motion of trapped vortices can be suppressed 
by making the superconductor films sufficiently narrow, 
it is currently not known how to suppress and control 
charge noise. Indeed, moving BCs of unknown origin are 
clearly responsible for the relatively short coherence time 
in the smaller Josephson charge qubits Q. In addition, 
BCs limit the performances of the flux qubits [3, the 
phase qubits [f| and the hybrid charge-flux qubits i0| be- 
cause they induce fluctuations in the critical current of 
the Josephson junctions and lead to qubit dephasing. 

It is firmly established now that in all superconducting 
qubits both the spectrum of charge noise and critical cur- 
rent fluctuations display a 1/f behavior at low frequen- 
cies. Specifically, charge noise spectra S q (f) = y- have 
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been directly measured up to frequencies / < 1 3 Hz 
00 0000 The intensity of the noise varies in 
the small range a = 10~ 3 — 10~ 4 e and because its mag- 
nitude is greatly reduced by the echo techniques [HJ, [Hi , 
it is commonly believed that 1/f noise does not extend 
above f m ~ 1 MHz. Similarly, it has been shown that 
critical current fluctuations display a 1// behavior at 
/ < 10 Hz alth oug h very little is known at higher fre- 
quency mm 13. 

Both charge noise and critical current fluctuations can 
be phenomenologically explained by modeling the envi- 
ronment as a collection of discrete bistable fluctuators, 
representing charged impurities hopping between differ- 
ent locations in the substrate or in the tunnel barrier. 
Because the hopping time is exponential in the height of 
the barrier separating different impurity states, one ex- 
pects a wide distribution of the relaxation times for the 
fluctuators. Such wide distributions naturally lead to the 
low frequency 1/f dependence. 

The discrete nature of the environment has important 
consequences for the qubit dephasing. As was shown in 
Ref. fl5L the effect of a bath of fluctuators on the qubit 
is very different from that of a continuum of linearly 
coupled oscillator modes, which is the subject of the 
Caldeira-Leggett dissipation theory 0. The discrete- 
ness of the environment becomes especially important 
for the slow modes which switching rates that are longer 
or comparable to the operational time of the qubits. 
Qubit interaction with slow fluctuators might lead to 
non-Markovian errors that are difficult to correct [TtI ] and 
to a fast initial qubit dephasing 0] ■ 

It has been recently suggested that dynamical decou- 
pling techniques |l9j| could be very promising strategies 
for suppressing the noise generated by slow fluctuators in 
superconducting qubits [2(3, l2ll I22L l23j . Unfortunately, 
because the behavior of the charge noise power spectrum 
is known only in very limited interval of frequencies, it is 
still not obvious whether it is possible to develop efficient 
dynamical decoupling schemes for these devices. 
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Indeed, besides the 1// low frequency part of the 
charge noise spectrum, so far it is only known that charge 
noise displays an ohmic behavior at very high frequen- 
cies (7 GHz < / < 100 GHz) The noise spectrum 
at intermediate frequencies can only be inferred from its 
effects on the qubit and so its detailed behavior is not 
known. This missing part is very important for the de- 
velopment of dynamical decoupling techniques, because 
the presence of noise at the time scales corresponding 
to the time interval between echo pulses limits the error 
suppression capability of these techniques. 

Interestingly, it has been argued that dynamical de- 
coupling techniques might be also used to extract infor- 
mation on this missing part of the noise spectrum, as 
a sort of diagnostic tools to infer spectral proprieties of 
the charge fluctuators coupled to the aubit|2a.l23|. Sev- 
eral numerical results pointed out that the interplay be- 
tween the applied pulses and the qubit dynamics might 
originate a rich behavior in the dissipative dynamics de- 
pending on the characteristics of the fluctuators and the 
working point of the qubit. This conjecture is supported 
also by a few analytical results. For example, in Ref. I22I 
the authors derived the analytical solution for the ideal 
(i.e. with no errors in the pulses) dynamical decoupling 
of a qubit coupled to many fluctuators at pure dephas- 
ing, while in Ref. the analytical solution for the ideal 
decoupling of a qubit coupled to a fluctuator at a generic 
working point has been sketched. However, we believe 
that from these analytic results it remains still difficult 
to estimate the effective diagnostic capabilities of the dy- 
namical decoupling techniques. Moreover, in a realistic 
set-up, one should expect that the dynamical decoupling 
is not ideal and it should be desirable to extend this anal- 
ysis to include errors in the pulses. 

This paper aims (i) to simplify the previous analytical 
solutions of the dynamical decoupling at pure dephasing, 
(ii) to find the analytic solution of the dissipative qubit 
dynamics at the optimal working point, (iii) to extend 
these analytical results to the case of non-ideal control 
pulses. We hope that this novel approach and our an- 
alytical solutions might be relevant in the design and 
interpretation of new experiments. 

A remark is needed at this point. In this paper we shall 
study only the dissipative dynamics of the qubit coupled 
to a single classical fluctuator that is subject to dynami- 
cal decoupling. There are several reasons for this choice. 
First, there are interesting experimental situations where 
there is evidence of a single fluctuator dominating the 
qubit dynamics 0, |2f|, |2(j . In this case, the analytical 
solutions derived in the following sections are sufficient to 
interpret the experimental data. Second, we believe that 
even from this simple solvable model one can obtain some 
important insights and intuitions on the slow frequency 
noise and especially on its non-Gaussian nature. Third, 
at pure dephasing, it is trivial to extend our analytical 
results to many fluctuators interacting with the qubit, 
because the total response is simply the linear combina- 
tion of the contributions from each individual fluctuator 



On the other hand, the fact that we are considering a 
classical fluctuator implies that our discussion most likely 
will be relevant to superconducting qubits made out of a 
Cooper pair box or more generally to those superconduct- 
ing circuits containing ultrasmall junctions (having area 
< 0.05/im 2 ) where the presence of few dominant classical 
fluctuators have been clearly observed in critical current 
fluctuations [27L . 

However, it is relevant to mention here that also exper- 
iments in phase qubits 0, l29j j and very recently in flux 
qubits |30j show the presence of single charged impurities 
strongly coupled to the qubit and seriously affecting its 
dynamics. It would be very interesting to sec whether 
dynamical decoupling sequences could help the charac- 
terization of the noise even in these cases. However, since 
it is reasonable to expect that those charged impurities 
are quantum Two Level Systems present in the insulat- 
ing barrier of the Josephson junctions |3lj . the analysis 
discussed here should be extended in order to take into 
account the quantum nature of the fluctuator. This will 
be the subject of future work. 

This paper is organized as follows. In Sec. HT1 we shall 
briefly recall the physics of the Cooper pair box, we shall 
derive the total Hamiltonian of the qubit in interaction 
with a bath of fluctuators and with applied an external 
control field and we shall finally explain how sequences 
of pulses can be generated in this device. In Sec. IIIII we 
shall present analytical results for the dynamics of the 
qubit subjected to ideal dynamical decoupling (i.e. with 
no errors in the pulses) at different qubit working points. 
Specifically, in Subsec. IIII Al we shall consider the qubit 
at pure dephasing and we shall greatly simplify the ana- 
lytical result given in Ref. I22I thereby immediately elu- 
cidating the efficacy of the dynamical decoupling as a 
diagnostic tool. In Subsec. IIIIBI we shall provide the an- 
alytical solution for the dissipative dynamics of the qubit 
subjected to dynamical decoupling at the optimal point. 
This result will allow us to comment on the evidence of 
the anti-Zeno effect in this regime [23, and to explain 
how the latter can be conveniently avoided. In Sec. I1VI 
we shall discuss dynamical decoupling sequences with im- 
perfect pulses. The consequences of imperfect bang bang 
pulses have been already discussed, for the case of a sin- 
gle fluctuator, in Ref. l32t Here, by using our approach, 
we are able to derive bounds on the pulse noise that can 
be tolerated by the device while still achieving dynamical 
decoupling. Finally, Conclusions are discussed in Sec. [V] 

II. COOPER PAIR BOX, CHARGE NOISE AND 
CONTROL 

A Cooper pair box Q, is a small superconducting 
island connected to a reservoir via two Josephson 
junctions (see Fig. The box is characterized 

by two energies: the Josephson coupling energy 

Ej ($ext) = \J E 2 n + E 2 ]2 + 2Ej 1 Ej 2 cos (ir$ ex t/$o) 
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FIG. 1: Sketch of a Cooper pair box afflicted by charge noise 
present in the amorphous materials of the substrate and the 
junction barriers). 

and the charging energy E c = e 2 /2C. Eji, i = 1,2 are 
the Josephson energies and <& e xt, $0 are respectively 
the external magnetic flux piercing the superconducting 
loop and the flux quantum. C is the total capacitance 
of the box and e denotes the electron charge. 

Under the assumptions that (i) the charging energy 
E c is much larger than the Josephson coupling Ej, (ii) 
the superconducting gap A of the box is larger than its 
charging energy (so to avoid quasi-particle excitations), 
(iii) the temperature T is small compared to all these en- 
ergy scales, it is possible to describe the dynamics of the 
Cooper pair box as a two state system, whose Hamilto- 
nian in the charge basis {|0), |1)} reads: 

SE Ej(* ext ) 
H s = — v.z H ^ Ux W 

where u x ^ z are the Pauli matrices. Both the bias 
5E C = E c (l — CgVg/e) and the Josephson energy can be 
easily tuned by varying, respectively, the applied gate 
voltage V g and the magnetic field & ex t threading the loop 
containing the island. 

The Hamiltonian given in Eq. JIJ describes the dy- 
namics of an ideal Josephson charge qubit. However, it 
is known that many charged impurities are present in the 
amorphous materials both of the substrate (on which the 
Cooper pair box is grown) and of the Josephson barri- 
ers. These charges, by interacting with the qubit, are 
responsible for its decoherence. 

In order to write an Hamiltonian for a realistic (with 
charge noise) Cooper pair box, one usually assumes that, 
irrespective of the exact microscopic nature, the charge 
impurities are charges that tunnel between two positions. 
Each charge is described as a quantum Two Level System 
(TLS) that is controlled by the Hamiltonian [H HI H3 

H T ls = er 2 + tr x + H env . (2) 



Here e is the energy difference between the two tunnel- 
ing minima, t is the tunneling amplitude and t z ^ x are 
the Pauli matrices acting on the impurity states. Each 
TLS is subject to dissipation due to its own bath with 
Hamiltonian H env . In this picture, the interaction be- 
tween several charged impurities and the charge states of 
the Cooper pair box can be written as: 

H S b = S £^v i T % z o z . (3) 

i 

By adding Eqs. QJ, J5J) and we finally obtain the 
Hamiltonian of a non-ideal (with charge noise) Josephson 
charge qubit. 

Fortunately, experiments performed on Cooper pair 
boxes indicate that this rather complicated model of the 
charge noise can be simplified. In fact, it has been ob- 
served that charge impurities that are responsible for the 
low frequency noise behave essentially classical, i.e. they 
become equivalent to classical fiuctuators each charac- 
terized by a switching rate and the coupling to the 
qubit Vi . In this approximation the interacting Hamilto- 
nian given in Eq. J3J simplifies to the following qubit- 
fluctuators Hamiltonian: 

H$ B = Y,vMt)a z (4) 

i 

and the quantum TLS dynamics given in Eq. Q is sub- 
stituted by the classical noise representing a sym- 
metric Random Telegraph Process (RTP) switching be- 
tween values ±1 with rate 7^ [3r3 |. If we assume a dis- 
tribution P(7) oc I/7 at 7 < 7a/ for the switching rates, 
we find for the total noise = J^i v i£i(t) acting on the 
qubit a 1/f spectrum S^l^l) = ^4/|w|, to < 7m, with co- 
efficient A proportional to the number of fiuctuators per 
noise decade rid and their typical coupling to the qubit 
(v) 2 jUI^. Note that we assume the switchings to be 
symmetric, that is, we assume the switching rate from 
the upper to the lower and from the lower to the upper 
TLS state to be the same. This is valid if the temper- 
ature is larger than the TSL level spacing so that there 
is no population difference between the two states. This 
restriction can be relaxed |l "J| but for simplicity we shall 
work in the symmetric (high temperature) limit. 

Therefore, in this classical picture, the realistic Hamil- 
tonian of a Cooper pair box afflicted by 1// charge noise, 
is given by the sum of Eq. and Eq. In the qubit 
eigenstates, we find that: 

AE x(^) 
H = —Pz + ~y~ [cos9p z - sinflpx] , (5) 

where AE = y/SE 2 + E'j, 9 = axcta.n(Ej/6E c ) defines 
the qubit working point and p x%z are the Pauli matrices 
in the rotated basis. 

The state of the qubit given in Eq. (JSJ) can be inter- 
preted as the state of a pseudospin in a static magnetic 
field B = (0, 0, AE) along z with noise both in the x and 
z directions. Manipulations of the Cooper pair box state 
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can be then achieved by applying an electromagnetic field 
which rotates in the x — y plane at frequency ow, at or 
near the pseudospin precession frequency AE . Given 
the Hamiltonian of the control radio-frequency (RF) field 
as follows: 

H r f(t) = fl(t) [cos(uj r ft + Lp)p x + sin(u} r ft + ip)p y ] (6) 

where <p is the phase of the RF field and f2(i) is the 
amplitude, the total Hamiltonian H tot of the Josephson 
charge qubit with charge noise and applied control field 
reads (in the rotating frame): 

H tot = AE ~ ur} ^ + Q ^ ^ cQs ^ + gin ^ 

+ [p z cos 8 + (p x cos(u! r ft) — p y sin(o; r /t)) sin 6\ 

In superconducting qubits typical values for the energy 
splitting are AE ~ 5 — 10 Ghz while the amplitude of the 
applied RF field Sl r f might extend up to 250 MHz. 

Different pulses sequences can be then generated by 
modulating the amplitude £l(t) in Eq. 0. Specifically, 
usual dynamical decoupling sequences are realized 
first by setting: 

N 

= n rf [&{t - jfer) - <&(t -kr- t p )] ; 
fe=i 

where N is the number of pulses, Tp is the length of each 
pulse, denotes the Hcaviside step function and r is 
the time interval between each pulse. Then, by choos- 
ing different values for the constant A = 2fl r fTp and the 
phase <p. For example, rotations by tt around the x or y 
axes, tt x and ity pulses, are achieved by fixing A = tt and 
by setting respectively (p = and ip = ir/2. Notice that, 
the so called bang-bang control, characterized by pulses 
that are infinitely short and strong so that there is no 
qubit dynamics during the time of the pulse, is achieved 
in the limit Tp — ► for fixed value of the constant A. In 
this paper we always consider only this limiting case. 



III. IDEAL DYNAMICAL DECOUPLING: 
EXACT SOLUTIONS 

In this section we solve analytically the dynamics of a 
Cooper pair box coupled to a single fluctuator and sub- 
ject to ideal (i.e. without errors in the pulses) dynamical 
decoupling both at pure dephasing and at the optimal 
working point for the qubit. 

The analytical solutions are derived by following the 
reasoning given in Ref. l40i The starting point is the 
Hamiltonian given in Eq. JSJ that describes a qubit in 
interaction with a single fluctuator. The charge noise 
is then given by \{t) = v£(t), where v is the coupling 
strength of the fluctuator and denotes a symmetric 
RTP characterized by the switching rate 7. 



It may be useful to recall here that in general, depend- 
ing on the ratio between v and 7, we can distinguish 
between two different types of fluctuators: those charac- 
terized by v < 7 (named weakly coupled or Gaussian fluc- 
tuators), and those with v > 7 (named strongly coupled 
fluctuators) that are responsible for the non-Gaussian 
features of the noise |4l|. The effect on the qubit dy- 
namics of weak and strong fluctuators is dramatically 
different. In fact, while the weak fluctuators are fast and 
produce essentially the same effects on the qubit as a 
bath of harmonic oscillators, thereby leading to the ex- 
ponential decay of the coherent qubit oscillations at long 
times, the strongly coupled fluctuators, being slow, are 
responsible for memory effects in the qubit dynamics, for 
a fast initial decoherence of the qubit signal and for the 
inhomogeneous broadening of the qubit spectrum. 

Due to the time dependence of the noise generated by 
the fluctuator, the Cooper pair box behaves as a pseu- 
dospin in a variable magnetic field. Therefore, solving its 
dissipative dynamics is a nontrivial task. Fortunately, be- 
cause of the bistable nature of the fluctuator, the problem 
can be simplified. In fact, at different times, depending 
on the state of the fluctuator, the qubit vector state pro- 
cesses on the Bloch sphere around one of the two effective 
static magnetic fields: B± — (±u sin 9, 0, AE ± v cos 8). 
As a result, a convenient strategy to study the decoher- 
ence of the qubit is following the evolution of the state 
vector on the Bloch sphere. In this picture, solving the 
dissipative dynamics of the qubit reduce to calculating 
the probability p(x, t) that at time t the state vector 
reaches a point x = [x, y, z) on the sphere. 

A. Dynamical decoupling at pure dephasing. 

At pure dephasing, 8 = 0, the analytic solutions for 
the dephasing rate T2 of a qubit interacting with a single 
fluctuator are well known both in the case of free induc- 
tion decay pi liol lilj and when sequences of tt pulses 
are applied to the qubit j23, • Here we derive these 
results by using a different approach that allows us (i) 
to simplify greatly the previous analytical solutions, (ii) 
to elucidate better the diagnostic capability of the dy- 
namical decoupling and finally (iii) to suggest an easy 
experiment that might be performed in order to extract 
information on the microscopic nature of the fluctuator. 

At pure dephasing the qubit acquires an extra random 
phase shift <j>(t) = v / * dt'x(t') due to the interaction with 
the fluctuator. The dephasing rate T 2 is then given by 
the decay rate of the average value: 

(e«)= J dcf>p(4>,t)e^ ■ (8) 

In order to evaluate it, we need to calculate the proba- 
bility distribution p((f>, t). 

Let us first introduce the small time steps ts = t/M 
and discretize the time integral given in Eq. JHJ . The ran- 
dom phase shift can be written as (j>{t) = vts J^ n=1 Xm 
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where Xn = x( n ^s) an( l the integration over time can be 
thought of as random walk process, where at each time 
step the random walker moves a step xs = vtg in a direc- 
tion depending on the current position of the fluctuator. 

To first order in t$, the probability that the fluctu- 
ator does not switch, i.e. that the next step is made 
in the same direction as the previous one, is given by 
a = Po(ts) = 1 — 7*s- Similarly, fi = Pi{ts) = yts is the 
probability that one switch occurs, i.e. that the step 
is made in the opposite direction of the previous one. 
Higher numbers of switches are negligible to the first or- 
der. 

Let us denote by m the number of steps from the ori- 
gin (i.e. the position of the random walker is x = mxs) 
and by n the number of temporal steps (i.e. the dimen- 
sional time is t n = nts). Let us calculate the discretized 
probability P n (m) that the random walker reaches the 
position to at time step n. It is useful to introduce the 
probabilities P^(m) and P~(m) for the random walker 
to reach respectively the position to coming from the left 
and from the right. These probabilities satisfy the fol- 
lowing master equations: 



P+ +1 (to) = aP+(m - 1) = pP~(m - 1) 
Pn+M) = PPn(m + 1) = aP-(m + 1) 



(9) 



and are such that P„(to) = P r f(m) + P~(m). 

In the limit of infinitesimal steps t$, these probabilities 
reduce to: 

p n( m ) = p{mx s ,nt s ) = P+(<t>,t) 
p n( m ) = p(mxs,nts) =p-(<f>,t) . 

Accordingly, the master equations given in Eqs. (j^J can 
be written as: 



dtp(<f>,t) = M deph p{(j),t) 
where we have introduced the matrix: 

-7 - vd<j, 7 



Mdeph = 



7 



and we have defined: p(</>, t) = 



-7 + vd^ 
p+(<M) 



(10) 



(11) 



M) /' 



The solutions of Eqs. (fTUf> can be easily guessed as 
p = Ce 4K< ^~ wt and we get the following dispersion equa- 
tion: 



where the matrix 



has the following eigenvalues 



-iuC = M deph C 



7 — IVK 7 

7 —7 + IVK 



w± = -i-f ± i-fJ 1 - -^k 2 = -£7 ± ijfj, n 



(12) 



(13) 



(14) 



and corresponding normalized eigenvectors: 
1 ( -i%K±fj, K 



1 



(15) 



It follows that the general solution to Eq. I|1U|I reads: 

r°° Hk 

p(<M) = / — [A + v+e- M +* + A_v_e- M -'l e iK4 
J -00 27r 



(16) 



where the constants A± are determined by the initial 
conditions: 



p(0) = A + v + +A_v_ 



(17) 



or, by using a more compact notation, A± — p(0) ■ v±. 
By introducing the matrices: 

U=( e -l.t) and y = (v + |v_) (18) 

we can conveniently write Eq. (|16|) as follows: 

r°° rih- 

V{ct>,t) = j -[A(i)p(0)]e^ (19) 

where k(t) = VUV T . Then, by substituting Eq. O 
into Eq. (JSJ) , we finally find that the average phase shift 
reads: 



(e**) = / dn I jp [A(t)p(0)]e 4 ( K+1 ^ 



d K [A(i)p(0)] 5{k + 1) 
[A(t)p(0)] K= _ x . 



(20) 



In order to evaluate the dephasing rate T2, it is conve- 
nient to expand the vector p(0) in terms of the eigenvec- 
tors of A(i). It is then a matter of simple algebra to find 
that the resulting decay of the average phase shift given 
in Eq. I|20|l is indeed given by the sum of two exponen- 
tials with coefficients = — In |A±|/i, where 

A± = e -^±W (21) 

are the eigenvalues of the matrix A(f) and we define 



1 _ m 

1 -.2 



The long time decay is given by the component that 
decays slower, i.e. the one that corresponds to the larger 
eigenvalue 

T 2 =-ln|max(A + ,A_)|/t = 7(l-3fyt) (22) 

This result agrees with the one given in Refs. mEEt 
A similar approach can be used in order to evaluate 
the dephasing rate I2 of the qubit that is subject to 
dynamical decoupling. In particular, let us assume that 
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we apply to the qubit a sequence of tt x pulses. It is easy 
to realize that the effect of each tt x pulse is simply to 
interchange the constants A + and As a result, when 
a sequence of N n x pulses is applied, the probability given 
in Eq. l(T§|l reduces to: 



rf/r 2 



P(0,t) 



2tt 



[A dd (t)p(0)] e 



(23) 



where: 




8 10 



A dd (i) =a x A(r)---a x A(r) 



(24) 



and t = Nt is the total time. 

Once this probability is known, it is straightforward 
to evaluate first the average phase shift acquired by the 
qubit: 



f dnf g [A dd (i)p(0)] 



dn [A dd (i)p(0)] S(k 



FIG. 2: Dephasing rate of a qubit interacting with a weak 
fluctuator as a function of the time interval r of the bang 
bang pulses. The analytical solution given in Eq. 128H is 
represented by the solid line, while the black dots indicate the 
numerical results. The values for the fluctuator parameters 
are respectively v = 0.1 and 7 = 1 (in unit of AE). The 
(25)dephasing ra te F 2 = v 2 /2^y is given in Eq. (J2HJ. 



[A dd (t)p(0)]_ 



and then the corresponding dephasing rate: 



■ln|max(A^,A^)|/i 



(26) 



Here X± are the eigenvalues of a x K(r). Depending on 
the nature of the fluctuator (i.e. weak or strong), the 
eigenvalues A± can be written as a function both of the 
fluctuator and the control parameters. 

Let us first consider that the fluctuator coupled to the 
qubit is weak, i.e. v < 7. We find that: 



1 and analytical results are in good agreement. Moreover, 
as one should expect, by reducing the time interval be- 
tween the pulses, the dephasing rate r dd vanishes. 

Let us now consider the case of a strongly coupled fluc- 
tuator. For v > 7, the constant /1 becomes complex 
(/i = ijl) and we find that the eigenvalues read now: 



sinli7/iT ±4/ cosh 7/iT ^ 

V 7 



(27) 



Since the eigenvalue A+ is the maximum one for each 
choice of the parameters, we find that the dephasing rate 
is given by: 



r dd 

1 2 



ln|A"|/t 



In I A, 



(28) 



In particular, to the lowest order in w/7 we find that: 



r 



dd 



V 

2T 



1 tan!i7T 

7T 



(29) 



From this simple analytical result, it is evident that in 
order to achieve an efficient suppression of the noise gen- 
erated by a weak fluctuator we need to choose the time 
interval between two pulses r < I/7. Namely, the ef- 
ficiency of the dynamical decoupling is related to the 
switching rate 7 of the weak fluctuator. 

In Fig. |2we plot r dd (r) given in Eq. l(29"|l and we com- 
pare our analytic solution with the result of a numerical 
simulation where the time evolution of the qubit was fol- 
lowed for a given realization of the fluctuator RTP and 
averages taken over many such realizations. Numerical 



A+ = 



sin7/ir ± \ / — — cos 2 7/ir 
V T 



(30) 



Differently from the expressions given in Eqs. I|27|l . here 
the two eigenvalues A± are oscillating in the time inter- 
val t between the pulses. It follows that there are two 
relevant dephasing times now: 



r dd± = -ln|A^|A. 
To the lowest order in 7/w, we find that: 



^dd± 



-ln|A±|/r = 7 



1 

1 ± — sin vt 

VT 



(31) 



(32) 



Notice that, by looking at Eq. (J33J), it is quite straight- 
forward to realize that for the case of strong fluctuator 
an efficient suppression of the noise is achieved by set- 
ting t < 1/v. Therefore, the efficiency of the dynamical 
decoupling is related to the coupling strength v of the 
strong fluctuator. 

In Fig.|3|we plot the dephasing rates r dd± (r) given in 
Eq. I|3I|) and we compare them with numerics. As evi- 
denced by Fig. |3 the results of the numerical simulations 
clearly show that, depending on the time interval r be- 
tween the pulses, the dephasing rate T dd of the qubit is 
dominated differently by the two dephasing rates T dd± . 

In order to understand this behavior it is useful to write 
the average phase shift given in Eq. I|25|) as follows: 



[N 



(33) 



7 



rr/r a 




FIG. 3: Dephasing rates of a qubit interacting with a strong 
fluctuate* as a function of the time interval r of the bang bang 
pulses. Solid lines represent the analytical solutions given in 
Eqs. (13111 while black dots are the numerical results. The 
values for the fluctuator parameters are respectively v = 0.1 
and 7 = 0.01 (in units of AE). From Eq. 1321 we have 
T2 = 7. 

where we have introduced the weights w±. These can 
be calculated as follows. First, we write the initial state 
p(0), given in Eq. I|17[) . as a linear superposition of the 
eigenvectors u± of the matrix a x K{r): 

i=± 

Then, we solve the following linear system of equations 
in order to find the coefficients Be 

p(0) • Uj = J2 Bi"* ■ fi J- ( 34 ) 

i 

For example, in the numerical simulations shown in 
Fig. |3 the fluctuator was initially assumed at equilib- 
rium, i.e. p(0) = (1, l)/2. By following the above proce- 
dure, we find that the weights read: 



2 2 v /cosh 2 1[ iT - £ 

In particular, to the lowest order in "f/v, we obtain that 
w±(t) = h(l ±cos?jt) and we can immediately explain 
the oscillating behavior displayed in Fig. [3] 

An important consideration is relevant at this point. 
By comparing Fig. [21 and Fig. [3] it is immediate to re- 
alize that the two plots are qualitatively different. One 
can then exploit this difference in order to distinguish 
between a weak and a strong fluctuator coupled to the 
qubit. For example, in an experiment, once a single flu- 
cutator has been singled out in a device, one could mea- 
sure systematically the dephasing time of the qubit after 
cycles of pulses characterized by different times r and 
observe if there is an oscillatory or not behavior in the 
data. The efficiency of the dynamical decoupling will 
then provide further information respectively on the cou- 
pling strength of the fluctuator v or its switching rate 7. 

However, in order to characterize completely the fluc- 
tuator both the two parameters v and 7 must be known. 



B. Dynamical decoupling at the optimal point. 

In this section we study the dissipative dynamics of a 
Cooper pair box tuned at the optimal point, i.e. 6 = 
that is subject to control sequences of it pulses. By using 
a similar approach to the one discussed in the previous 
section, we find the analytic solution for the qubit deco- 
herence rate. 

Contrary to the case of pure dephasing, at the opti- 
mal point the performance of the dynamical decoupling 
depends both on the characteristics of the charged fluctu- 
ators coupled to the qubit and on the qubit energy level 
splitting. In particular, quite worryingly, it has been no- 
ticed that the qubit decohcrence can be even accelerated 
if the cycle time of the control is long with respect to the 
period of the free evolution of the qubit ■ This 

acceleration is sometimes referred in the literature as the 
anti-Zeno effect [42^ . 

In the following we shall see indeed that this accelera- 
tion appears only when sequences of ir x pulses are used 
to flip the qubit state and can be conveniently avoided 
by flipping the state with w y pulses. 

To this aim, let us introduce the time evolution oper- 
ator U(t) = Te'^o H (*')<**', where H(t) is the time de- 
pendent Hamiltonian given in Eq. and let us study 
the time evolution of the qubit when a cycle of tt x and 
ir y pulses is applied to it. Let r be the time interval be- 
tween the pulses and let us consider first the case of pure 
dephasing. 

Since the time evolution operator U(t) depends only 
on the qubit operator a z , at pure dephasing we see that 
the noise coming from fluctuators that are static during 
the time interval r can be completely eliminated either 
by applying a cycle of n x or ir y pulses. In fact, mathe- 
matically we obtain that: 

a x U(T)a x U(T) = a y U(T)a y U(T) = 1 (36) 

where 1 is the identity. 

On the contrary, at the optimal point, the evolution 
operator U(t) depends both on the qubit operator a z 
(because of the interaction with the fluctuators) and a x 
(because of the Josephson coupling energy Ej). Mathe- 
matically we find that: 

(T x U(t)ct x U(t) ^ 1 and a y U (t)(T v U (t) = 1 (37) 

Namely, only by flipping the qubit state around the y-axis 
we can completely suppress the noise due to fluctuators 
that are static in the time interval r. Therefore, at the 
optimal point, the performance of the dynamical decou- 
pling changes remarkably depending whether sequences 
of 7T X or TT y pulses are used to decouple the qubit. 

These considerations are clearly illustrated in Fig. 0] 
and Fig. |SJ In Fig. 0] we study numerically the dissi- 
pative dynamics of a Cooper pair box that is coupled 
to a weak fluctuator. The qubit dephasing rate is 
displayed as a function of the time interval r between 
cycles of ir y pulses (upper Panel) and ir x pulses (lower 
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FIG. 4: Dissipative dynamics of the dynamical decoupling 
applied to a Cooper pair box coupled to a weak fluctuator at 
the optimal point (Ej = 10, v = 0.1, 7 = 1). From Eq. 1561 
we have T2 = ^v 2 /Ej. The lines are added to guide the eye 
and do not represent any analytical result 
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FIG. 5: Dissipative dynamics of the dynamical decoupling 
applied to a Cooper pair box coupled to a strong fluctuator 
at the optimal point (Ej = l,v — 0.1,7 = 0.01). From Eq. 
we have F 2 = -yv 2 /Ej. 



Panel). Notice that only when tt v pulses are applied to 
the qubit, the dephasing rate is monotonously re- 
duced as the nipping rate is increased. On the contrary, 
for sequences of ir x pulses, there arc values of r where the 
dephasing rate increases remarkably as compared to 
the rate when there are no control-flips (i.e. anti-Zcno 
effect HIH). 

A similar analysis is performed in Fig. [5] where the dis- 
sipative dynamics of the qubit coupled to a strong fluc- 
tuator is illustrated. Again, we find that the dephasing 
rate displays both the anti-Zeno effect and rapid os- 
cillations when sequences of tt x pulses are applied to the 
qubit (lower Panel). However, both the the acceleration 
and the oscillations disappear when control sequences of 
n y pulses are instead used (upper Panel). 

Note that when control sequences of ir x pulses arc used 
both in Fig. 0] and Fig. [5] the efficiency of the dynamical 
decoupling is related to the qubit energy scale Ej . In this 
case, in fact we need to choose time interval r < l/Ej in 
order to obtain an efficient suppression of the noise due 
to the fluctuator. For sequences of tt v pulses the effi- 
ciency of the dynamical decoupling is determined cither 
by the switching rate r < 1/7 (if 7 < Ej) or energy scale 
t < l/Ej (if 7 > Ej) for the weak fluctuator and by 



r < I/7) for the strong fluctuator. This will be discussed 
after we have derived the analytic expressions below. 

As a final remark, let us notice that the different per- 
formances of sequences of tt x and n y pulses can be easily 
explained in the context of effective Hamiltonian theory 
[43j, by stating that while ^-pulses provides decoupling 
only to the first order in the Magnus expansion of the 
dephasing rate, 7r y -pulses gives exact decoupling to all 
orders. 

Having clarified this point, let us now proceed to de- 
rive the analytical solution for the dissipative dynamics 
of the qubit at optimal point. In view of the above con- 
siderations, in the following we shall only consider control 
sequences of Tr y pulses. 

The analytical solution is derived by following the same 
reasoning discussed in the previous section. Here, the 
main difference is that the two effective static magnetic 
fields around which the state of the Cooper pair box pre- 
cesses read now: Bj- = (±v, 0, Ej). As a consequence, at 
the optimal point, the qubit does not acquire only a sim- 
ple random phase due to the interaction with the charged 
fluctuator, but its state will spread over the entire Bloch 
sphere. 

In order to calculate the probability p(x, t) that the 
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qubit vector-state reaches the point x = (x,y,z) on the 
Bloch sphere at a certain time t, first we split the proba- 
bility p(x, t) = p + (x, t) + p_ (x, t) and then we derive the 
master equations for p + (x, t) and p- (x, t) : 



p+(x, t + t s )= ap+{U- x x, i) + 0p_(U+ l x, t) 
p-(x,t + t s ) =^p + (UZ 1 x,t) +ap-(UI 1 ic,t) 



(38) 



where a = 1 — 7is is the probability for the fluctua- 
tor to remain in the same state at the time ts while 
= its denotes the probability that it switches. No- 
tice that in Eqs. (|38[) we have introduced the operator 
U± = e* sB±R , where R = (R x ,R y , R z ) and the rotation 
matrices are defined as follows: 




0-1] 
1 



1 
| 
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By performing the continuous limit ts — > 0, from Eqs. 
(|38|l we obtain directly the equations governing the time 
evolution of the Cooper pair box before the first n v pulse 
is applied: 



dtP+ = IP- - ( B + ' R- x ) • V x p+ - 7P+ 
dtP- = 1P+ - (B_ ■ Rx) • V x p_ - 7p_ 



(39) 



In this picture to apply a w y pulse to the qubit simply 
means to reverse the effective magnetic field. As a result, 
given Eqs. (|39|) . it is immediate to write the equations 
describing the time evolution of the state of the Cooper 
pair box after the application of the first ir y pulse: 



dtP+ = 1P- 
dtP- = 1P+ 



(B+ • Rx) • V x p+ - 7p+ 
(B_ • Rx) • V x p_ - 7P _ 



(40) 



Clearly, once the second ir y pulse is applied, the dissipa- 
tive dynamics of the box is again described by the set of 
equations given in Eqs. (|39[) and so on for all the control 
cycles. Still we need to understand how to match the 
solutions just before and after the application of each ir y 
pulse. 

Let p±(x, t — > 0~) be the solutions of Eqs. (|39f) just 
before the first pulse, while p±(x, r + ) are the so- 
lutions just after the pulse. Notice that the latter are 
indeed initial conditions for Eqs. (|40|l . It is then easy to 
realize that a 7r rotation around the y-axis leads to the 
following matching conditions: 



P+(x,r -> 0+) =p+(x,r -> ) 
p-(x,r 0+) =p_(x,r -> 0") 



(41) 



Therefore, we can write explicitly Eqs. (|40|l as follows 
d t p+ = [Ej(yd x - xd y ) + v(zd y - yd z )]p + - jSp 



d t p- = [Ej(yd x - xd y ) - v(zd y - yd x )]p- + jSp 
where dp — p+ — p~. 



(42) 



In order to solve them, it is useful to define the follow- 
ing expectation values: 



k± = J d 3 kp±(x 7 t)k 



x,y,z 



(43) 



We can then write Eqs. I)42|l in terms of Eqs. I|43|l and we 
obtain the following set of equations 



d t x + = -7 (x + - x_) - Ejy + 

dty+ = -7 (y+ -y-) + e .j%+ ~ vz + 

d t z + = -i(z+ - 2_) + vy + 
d t x- = 7(a; + -xJ)- Ejy_ 

d t y- = i(y+ -y-) + Ejx_ + vz_ 

d t z- = 7(2+ -z-)- vy- 



(44) 



having constant coefficients. Moreover, by introduc- 
ing the variables X± = x + ± x_, Y± = y+ ± y_ and 
Z± = z+ ± z-, Eqs. I|44|l can be further reduced into the 
following two independent sets of equations: 



8 t X+ = -EjY+ 
t Y + = EjX+ - vZ- 
dtZ- = -2 7 Z_ + vY+ 



(45) 



and 



d t X_ = -2 7 X_ 
d t Y_ = -2 7 y_ - 
d t Z+ = vY_ 



- EjY- 
EjX. 



vZj, 



(46) 



Finally, by solving the system of Eqs. I)45[l. we can esti- 
mate the dephasing rate Yi. To this aim, it is sufficient 
to prepare the system in an equal superposition of eigen- 
states, i.e. with Z = 0, and then to study the decay of 
the x and y components of the Bloch vector. Namely, we 
need to derive the decaying behavior of the variables X + 
and Y + . 

It is convenient to write the system of Eqs. (|45|l in 
matrix form: 



(47) 











9 A 


'y + 








\ Y z + -i 







where the matrix M\ reads: 



Mi 






-Ej 





Ej 





— V 





V 


-2 7 



(48) 



After a w y pulse is applied to the qubit, the matrix Mi 
is replaced by the following matrix: 



Mo 



Ej 
-Ej V 

-v -27 



(49) 
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By resorting to Eqs. Q41JI. it is straightforward to see 
that the evolution of the qubit state after one full control 
cycle can be described by the operator: 

T = n y e AU - T n y e MlT (50) 

where we introduced the matrix IT,, = diag(l, 1,1). 
Eq. (|50|l can be further simplified to: 

T=(U y Le M ^f (51) 

where we defined the matrix L = diag(l, —1,1). Finally 
we find that, after an echo cycle: 

3 

Wl e- T i * i =1,2,3 (52) 

i=l 

I 



where 

if^-lnH/r (53) 

and fj,i are the eigenvalues of the matrix T. It follows 
that the dissipative dynamics of the qubit tuned at the 
optimal point is controlled by three different dephasing 
rates. 

Let us now discuss the decay rates in the limiting cases 
of weak (v <C 7) and strong (v 7) fluctuator. In 
both cases, it can be proved that in Eq. (|52|l the weight 

u'3 =0. As a result, only the two dephasing rates 

and r^ 2 ^ are relevant for the qubit dissipative dynamics. 

Let us first evaluate them for the case of a weak fluc- 
tuator. After some calculations we find that: 



-dd(l) 



Ej + 4 7 2 



Ej - 4 7 2 sin (Ejt) 



Ej 



47 2 Ejt 



8 7 2 



Ej + 4 7 2 



Ejt\ tanh (77") 



(54) 



ndd(2) 



Ej + 4 7 2 



El 
Ej 



47 2 sin (Ejt) 
47 2 Ejt 



8 7 2 



Ej 



4 7 2 



Ejt\ coth (77-) 



1 T 



(55) 



Both from Eq. lf5^|l and Eq. ljSl)|) it is evident that, 
contrary to the case of pure dephasing, the energy level 
splitting of the qubit Ej plays now a major role in the 
dissipative dynamics of the qubit. 

In order to understand the efficiency of the dynamical 
decoupling, it is useful to discuss the dephasing rates in 
two limiting case. 

(i) weak fluctuator with 7 < Ej . 
In this case Eq. (|S"4fl and Eq. lf53|) reduce to: 



larger weight is given to this solution. 



rf/r 2 



r 



dd(i) 



r 



dd(2) 
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1 - 


sin (Ejt) 




Ej 


Ejt 




2 


1 + 


sin (Ejt) 


8 7 2 


Ej 


Ejt 


Ej 



(56) 
coth (77-) 



1 T 




From Eq. 1|56[) it is evident that in order to achieve an ef- 
ficient suppression of the noise induced by the fluctuator 
we need to apply sequences of n y pulses with time inter- 
val between the pulses r such that t < ir/Ej. Once this 

constrain is satisfied, we find that T 



dd(2) 



71?' 



2 ( 7 t) 2 . 



2 r ~ B 2 j 3 

Since this represents a negligible term, the decay of the 
Cooper pair box is effectively dominated by the decay 
rate rf (1) . 

As an illustrative example, in Fig. we consider the 
dissipative dynamics of a qubit that is coupled to a fluc- 
tuator in this regime and we compare the exact solutions 
given in Eq. (|54f) and Eq. (|55|) to the numerics. Notice 
that for small time interval r the dissipative dynamics is 
governed by the larger decay rate, presumably because 



FIG. 6: Weak fluctuator (Ej = 1, v = 0.1, 7 = 1). The 
lines show the Eq. 11541 and Eq. 1551 . points are numerical 
simulations. From Eq. I54H we have V2 = yv 2 /(Ej + 47 2 ). 

(ii) weak fluctuator with 7 3> Ej. 
In this limit, Eq. ijSH and Eq. (|55)l become: 



^dd(l) 



^dd(2) 



V 

4- 

4- 



sin (Ejt) 2 ( Ejt\ tanh (77-) 

_L I w _■ COS 



1 - 



EjT 

sin (Ejt) 

EjT 



1 T 



(57) 



The efficiency of the dynamical decoupling can be easily 
inferred from the decaying rate r^ 2 "* given in Eq. I|57|) 
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and, as in the previous case, the condition to satisfy is: 
t < n/Ej. By expanding in the limit Ejt <lwe obtain 

that rf (1) a £ (2 - - l) . It follows 

that in order to keep the latter negligible one has also to 
satisfy the condition r « 1/7. As a result, in this limit, 



one should choose as condition for the efficient decou- 
pling the strongest one, i.e. t <C I/7. This behavior is 
illustrated in Fig0 

Let us now consider the case of a strong fluctuator 
coupled to the qubit. The dephasing rates read now: 



-,dd(l) 



A 2 " 

IE] 



sin (rA) 



2A 2 



-A 



2 cos(rA) + 1-3 



sin(rA) 



A 2 + Ej + v 2 cos(rA) 



4^ [cot(^)-^]- 
E 2 +A 2 cot 2 (^) 



(58) 



2£ 2 



1 



, sin(r A) 
' rA 



A 2 + Ej + v 2 cos(rA) 



where A = y/E 2 , + v 2 . Since we are only considering 
noise that is weak compared to the level splitting of the 
qubit, v <C Ej, and the fluctuator is strong, i.e. 7 < v, 
only the case 7 <C Ej is relevant. In order to have effi- 
cient dynamical decoupling of ' it is immediate to 
realize that we need to set r < n/Ej. Once this con- 
strained is satisfied, we find that T2 d(2) « 0. Therefore 

the dissipative dynamics is finally dominated by . 
An example is shown in Fig. [SJ 



rf/r 2 




FIG. 7: Weak fluctuator (Ej = 1, v = 0.1, 7 = 100). The 
lines illustrate Eq. I|54^ and Eq. 155^ and the black points 
represent the numerical results. The inset shows the analytic 
and numerical solutions for small values of r. From Eq. 1571 
we have F 2 = v 2 /Ary. 



IV. IMPERFECT PULSES. 

So far we have assumed that the control pulses are 
ideal. Namely, they are zero-width 7r pulses applied ex- 
actly along the x or y axis of the Boch sphere. However, 
in real applications it is plausible that the duration of 
the pulses exhibits random fluctuations, leading to error 




2.5 5 7.5 10 12.5 15 17.5 20 

FIG. 8: Strong fluctuator (Ej = 1, v = 0.1, 7 = 0.01). The 
solid lines show Eq. 15411 and Eq. (I55t and the black points 
illustrate numerics. From Eq. 1581 we have F 2 = 7« 2 /Ej. 



in the rotation angle and that the effective rotation axis 
might be deviated from the x or y axis. 

Clearly, if the pulses are imperfect and the repetition 
rate between the control cycles is high, there is the pos- 
sibility that the dynamical decoupling introduces more 
dephasing than what it is effectively capable to remove. 
As a result, given the noise characteristics of both the 
environment and the control pulses, there must exist an 
optimal pulse rate that sets the best decoupling perfor- 
mance. This section aims to determine this optimal pulse 
rate. 

Control sequences with imperfect pulses have been al- 
ready considered in two recent works |32j, |44 | . Our anal- 
ysis differs from the previous studies because: (i) it takes 
into account the qubit dynamics between the control 
pulses; (ii) by using the analytic solutions derived in the 
previous sections, we are able to provide bounds on the 
amount of noise that can be tolerated in the pulses while 
still achieving efficient decoupling. 

In order to calculate the optimal pulse rate it is conve- 
nient to proceed by steps. First we study the dissipative 
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dynamics for an ideal qubit (without charge noise) that 
is dynamically decoupled with control sequences of noisy 
pulses, then we extend this analysis to a qubit afflicted 
by a charge fluctuator. 

A. Dephasing due only to errors in the pulses 

In this section we study the dissipative dynamics of 
an ideal qubit that is subject to dynamical decoupling 
with noise in the n pulses. Specifically, we assume that 
the effective rotation realized by each pulse fluctuates 
around the expected values: i.e. 6 — ir ± 8, where 5 is a 
random variable. 

In the absence of a fluctuator the time evolution of 
the vector state on the Bloch sphere from the end of one 
pulse to the end of next pulse is then given by 

Ui = e (^+Si)R, e SE e rR z at pure dephasing 
Ui = e ^+ s *) R y e E -'T R * at optimal point 

After ./V pulses, we find that the total time evolution 
reads: U = YliLi Ui- Let us now assume that the fluc- 
tuations Si (i = 1, . . . N) are uncorrelated and that they 
are Gaussian distributed. The average time evolution 
can then be easily calculated as (U) = Yli(Ui) — (U) N 
and the decay rate of the qubit can be inferred by eval- 
uating the eigenvalues of (Ui). Both for the case of pure 
dephasing and at the optimal point we find that, in the 
limit (<5 2 ) <C 1, the eigenvalues are: 

= -Pc 



a± = - cos (At) (1 - p c ) ± y - cos 2 (At) (1 - p c ) 2 + p c 

where p c = (cos 8) and A = 8E C at pure dephasing and 
A = Ej at qubit optimal point. 

It follows that, to the lowest order in (8 2 ), the dissipa- 
tive dynamic of the qubit is controlled by three possible 
decaying rates rlj = — In «i|/r for i — 0, ±. Precisely 
we find that: 

r dd(0) = <g> (5Q) 

r dd(±) = (0 (1±cosAr) (60 ) 

At long times we expect that the qubit decay is controlled 
by the largest dephasing rate given in Eq. |J5!5J. 

At small times, the controlled dissipative dynamics is 
more complicated. In order to understand it, we con- 
sider the case of a qubit that is tuned at pure dephasing 
AE = Ej = 1 and it is subject to sequences of ix x pulses 
with fluctuating noise —0.12 < 8 < 0.12 uniformly dis- 
tributed. We solve numerically the dissipative dynamics 
for two different initial preparation of the qubit: either 
the qubit initial state is prepared along the x-axis or the 
y-axis of the Bloch sphere. The results of our simulations 
are illustrated in Fig. |5J 
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FIG. 9: The three decay rates ' (black), (blue), 
rig (green) (8 uniformly distributed between -0.12 and 0.12 
and A = 1). The black (blue) points are the result of numer- 
ical simulation, with the initial state on the y-axis (a-axis) 
of the Bloch sphere and pulses around the a;-axis. The inset 
shows the weights (black) and A- (blue) for the initial 
state on the y-axis. If the initial state is on the s-axis the 
weights interchange. 



Notice that, depending on the initial conditions, the 
dissipative dynamics is dominated by different dephas- 
ing rates. For a qubit initially prepared along the y-axis 
(x-axis) of the Bloch sphere the dynamics is dominated 
by the dephasing rate r 2 + (r^ - )- In order to under- 
stand this, we need to evaluate the weights of the different 
components. 

Let us indicate with So the initial qubit vector state. 
It is convenient to write it as a linear combination of 
the eigenvectors v, = (v®, vf , v^) of the average oper- 
ator (Ui) describing the evolution of the qubit vector 
state on the Bloch sphere during a control cycle, i.e. 
So = Y^ K ^ v ii with k = 0, +i — ■ After N control cycles, 
the qubit vector state reads: s = (U) N s = J2 K A K a^ 
and the coefficient Ai can be easily found. 

In the considered example, when the qubit is ini- 
tially prepared along the y-axis of the Bloch sphere, i.e. 
So = e y , we find that the coefficients reads: Aq = 
and A± = (1 ± cos At) /2. Therefore we understand im- 
mediately why in Fig. [5] the dephasing rate r^ ^ never 
dominates the dissipative dynamics. Moreover it is clear 
that, according to the oscillatory behavior of the coeffi- 
cient A± displayed in the inset, the dissipative qubit dy- 
namics at different time intervals r between the control 
pulses is alternatively dominated by the two dephasing 
rates r^ d± . 

Notice that, in a realistic situation where dynamical 
decoupling sequences are introduced in the protocol dur- 
ing the computation, we might not know exactly which 
is the effective state of the qubit. As a result, from the 
point of view of the efficiency of the control when fluc- 
tuating noise in the pulses is considered, it is reasonable 
to assume that the dissipation is always given by the 
largest decay rate r^*' '' at all interval times r between 
the pulses. 
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B. Dephasing due to error in the pulses and 
fiuctuator noise 

Let us now consider that during the dynamical decou- 
pling there is both noise coming from the pulses and a 
charge fiuctuator that is coupled to the Cooper pair box. 
Since the two sources of noise are independent, we ex- 
pect that the relaxation rates will just add. We could not 
prove this conjecture analytically, but results of numer- 
ical simulations clearly demonstrate that this is always 
the case. 

As an example, in Fig. HOI we illustrate the dissipative 
dynamics of a qubit that is coupled to a weak fiuctuator 
at pure dephasing. In the figure we plot: the analytic 
solution for the dissipative dynamics when there is noise 
only in the pulses given in Eq. l(5§|l (green), the analytic 
solution for the dephasing rate due to the interaction of a 
weak fiuctuator derived in Eq. (|29() (blue) and the total 
dephasing rate resulting from the sum of the two (black) . 
Notice that there is a good agreement between numerics 
and analytic results. 




T 



1 2 3 4 5 



FIG. 10: T2 with 5 uniformly distributed between -0.04 and 
0.04 and A = 1, Ej — 0. One weak fiuctuator, v = 0.1, 
7 = 1. Green - noise only in the pulses, blue - noise only from 
the fiuctuator, black - the sum of the two noise contributions. 
Black points are numerical simulation. From Eq. 129^ we have 

r 2 = v 2 /2~f. 



the order of the rate without dynamical decoupling) there 
will be a range of r below this point where dynamical 
decoupling is working, and it will not be destroyed by 
the pulse noise. 

In Sec. IIIII we have demonstrated that, depending on 
the nature of the fiuctuator, efficient dynamical decou- 
pling is achieved for r < I/7 if the fiuctuator is weakly 
coupled to the qubit and for t < 1/v if the fiuctuator 
is non-Gaussian. As a consequence, depending on the 
type of fiuctuator interacting with the qubit, we find that 
the control sequences are efficient when the noise in the 
pulses satisfy the following conditions: 




for weak fiuctuator 



(S 2 ) < — for strong fiuctuator 

v 

Notice that the noise tolerated in the pulses is bounded 
by the second power of the ratio v/j < 1 for a weak fiuc- 
tuator and only by the first power of j/v < 1 for a non 
Gaussian fiuctuator. Therefore, in a qubit afflicted by a 
strong fiuctuator, efficient dynamical decoupling can be 
achieved with higher noise in the pulses. 

A similar analysis can be done for a qubit that is tuned 
at the optimal point. Using the analytic solutions given 
in Sec. IIII Bl we find that an efficient dynamical decou- 
pling can be achieved 

- For weak fiuctuator, i.e. v < 7, when: 

(S 2 )<^ if 7 «£./; 
(5 2 )< ^ if 7 »£./- 

- For strong fiuctuator, i.e v > 7, when: 



V. CONCLUSIONS. 



By looking at Fig. EH we see that, by decreasing the 
time interval t between the pulses, the dephasing rate ini- 
tially decreases, then at a particular value of t it reaches 
a pronounced minimum where the best dynamical de- 
coupling performance is achieved and finally, by further 
decreasing r, the performance of the decoupling deteri- 
orates. Since the efficiency of the decoupling (i.e. the 
value of the minimum) is related to the level of noise in 
the pulses, it might be interesting to estimate how much 
noise in the pulses can be tolerated in order to have an 
efficient suppression of the noise due to the fiuctuator. 

To this aim, one has to find the time interval r between 
the pulses when dynamical decoupling starts to reduce 
the decoherence rate. If the decoherence rate because of 
the noise on the flips (~ (S 2 )/t) at this pulse rate is small 
compared to the one from fluctuators (which will be of 



In this paper we studied the dissipative dynamics of a 
Cooper pair box (qubit) afflicted by charge noise that is 
dynamically decoupled. By focusing on the case where a 
single fiuctuator is coupled to the qubit, we investigated 
the efficiency of ideal (without errors in the pulses) de- 
coupling sequences. At qubit pure dephasing, irrespec- 
tive of whether one applies sequences of ir x or n y pulses, 
we found that the efficiency of the decoupling is strongly 
dependent on the characteristics of the fiuctuator. For 
the case of a weak (Gaussian) fiuctuator an efficient sup- 
pression of the noise is achieved by choosing the time 
interval between the pulse less than the switching rate 
of the fiuctuator, while for a strong (non-Gaussian) fiuc- 
tuator the relevant decoupling rate is determined by the 
coupling strength of the fiuctuator to the qubit. As a 
result, in this limit dynamical decoupling techniques can 
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be successfully implemented as a diagnostic tool to in- 
fer spectral informations on the nature of the fluctuator 
coupled to the qubit. 

A different situation is encountered at the qubit op- 
timal point, where we showed that the decoupling per- 
formance varies considerably depending whether control 
sequences of tt x or ix y pulses are used to decouple the 
qubit. Both for the case of weak and strong fluctuator 
we found that the dissipative dynamics is dominated by 
the qubit energy scale. Although for the special case of 
weak fluctuator characterized by a switching rate greater 
than the qubit energy level splitting we observed that 
the efficiency of the dynamical decoupling is achieved by 
choosing time interval between the pulses less than the 
switching rate of the fluctuator. 

We have extended this analysis in order to include also 
errors in the pulses. We showed that, depending on the 
characteristic of the fluctuator and the error in the pulses, 
there must exist an optimal pulse rate that sets the best 
decoupling performance. Using the analytic solutions for 
the dissipative qubit dynamics derived both at pure de- 



phasing and at the optimal point, we found upper bounds 
on the level of noise present in the pulses that can be 
tolerated while still achieving an efficient dynamical de- 
coupling. 

We believe that an interesting future research direc- 
tion is to extend this analysis in order to take into ac- 
count the quantum nature of the fluctuator. It would be 
extremely interesting to see whether dynamical decou- 
pling techniques could be useful in order to understand 
the microscopic mechanisms of relaxations at low Tem- 
perature of the quantum Two Level Systems present in 
the amorphous materials of the barrier and substrate. 
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